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Abstract 

We discuss the leptoquark interpretation of the anomalous HERA positron-jet 
events in the context of the YFS exponentiated Monte Carlo event generator treat- 
ment of the attendant multiple photon radiative effects for both the would-be signal 
and the SM background, wherein finite-py photon effects are properly taken into 
account and wherein infrared singularities are cancelled to all orders in a. We show 
that the HI and ZEUS data are consistent with such an interpretation for a lep- 
toquark coupling < 0.3gw, mass ~ 200 GeV, and width < 2GeV. Possible future 
tests are proposed. 
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Recently, the HI and ZEUS Collaborations have reported |]l|, |^ an anomalous amount 
of high , high x Bj events in the deep inelastic e"'"p scattering at HERA. In this paper, 
we investigate in detail the interplay between the rigorous treatment of the attendant 
multiple photon radiative effects, which we treat by the YFS P exponentiated Monte 
Carlo technique introduced in Ref. @] and applied to deep-inelastic lepton-nucleon scat- 
tering in Refs. [|, ^, see also Ref. 0, in the event generator LESKO-YFS, and the HERA 
phenomenon if it is interpreted as the exchange of a spin 1 or object in the s-channel in 
the reduced e'^^q^ parton level process at HERA. For, to our knowledge so far, only struc- 
ture function based, zero-pr, treatments of the higher order QED radiative effects have 
been applied to the HERA data fl], ^ so that our work will actually serve two purposes: 
(1), it will show that the proper treatment of the higher order finite-pT photon effects in 
the multiple photon radiation cannot account for the phenomenon; and that, (2), when 
a resonance of the appropriate coupling and mass is introduced, the HI and ZEUS data, 
in the presence of the realistic ^(7) radiation do indeed agree with the theoretical ex- 
pectations. Of course, other possible interpretations of the data are also possible and we 
refer the reader to Refs. p, |, for further discussion of such approaches to the HERA 
phenomenon. 

For definiteness, we shall call this object a 'leptoquark', since according to Ref. 



it has a mass ~ 200 GeV and a width < 2.2 GeV. We stress as it has been already 



done in Ref. [12| that as long as the couplings of this object, which may be composite or 
elementary, are sufficiently chiral and diagonal in flavor and zero on diquark fields, it is 
not in contradiction with any known physical requirement. In Ref. [0, we have discussed 
the constraints implied by recent LEP2 data on this interpretation of the HERA data 
due to crossing it into the e~^e~ — » hadrons channel, as well as the constraints derived 



in Refs. [|T5|, |T6| . See also Refs. @, |T^, |l^ for related and similar analyses of the crossed 



reaction and see Refs. |TD[ for a discussion of contraints on this interpretation implied 

by recent TEVATRON data. The net result of these discussions is that the leptoquark 
interpretation is still viable if the coupling strength to fermions is small enough, < 0.3(7vy, 
where gy/ is the EW coupling constant. Thus, we proceed entirely phenomenologically 
and try to answer the very definite question as to whether the HERA data are consistent 
with a 'leptoquark' explanation if we treat the attendant higher order radiative effects 
via the YFS exponentiated Monte Carlo methods of two of us (S.J. and B.F.L.W.) as 
they are realized in the event generator LESKO-YFS p, ^, for both the signal and the 
background. 

More specifically, we first record the differential cross section for e^^q' — > e'^'^cp in the 
presence of the HERA leptoquark for two models of leptoquark charges. We have, for the 
outgoing e"*" elemental solid angle dVLe+ in the CM system, 

^ \M\\ (1) 



where the spin averaged squared matrix element, for the leptoquark (X) charges (Q, Q — 1) 
with Q = 5/3 coupling to left(right)-handed quarks {model (l)(model (2))} and with 
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Q = 4/3 coupling to left(right)-handed quarks {model (3)(model (4))}, is 
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where h{H) = — 1(+1) for the left [H = L)-, right [H = i?)-handed cases respec- 
tively, S is the spin of X, 6ab is the Kronecker delta function, uq = 5q -^u + 5q ±s, 
sq = 6q 5S + 6q iu and where we have defined the following kinematical and dynami- 
cal variables: Dx{s) = s - + iTxs/Mx, Mx = 200 GeV, Tx <2 GeV, Dz{i) = 
i-Ml DA{i)=i,vzif) = lh-QfSm^9w, az{f) = VA{f) = Q aM)=0, 
Gz = e/{sm9wcos9w), Ga = e, gx = (e/ smew){l + 5), s = (pi + qif, i = (pi - 
P2Y, u = {Pi — ^2)^5 where {^1,^1} are the incoming {e~^,^q^ 4-momenta respectively 
and {^2,92} are the outgoing {e"*", ^q^} 4-momenta respectively. Here, 6 is unknown and 
is to be varied to see what the data will allow. (According to Refs. [|13], [I^ 
crossed reaction constraints we hope to find that —1.5 ^ S ^ —0.5 for spin 1. 



!], from the 
< —0.5 for spin 1.) Is is the 
usual weak isospin 3-component for fermion / and Qf is its electric charge in units of the 
positron charge e. We have thus complied with the constraint from Refs. [|12|, O, |16[ that 



only quarks of a specific chirality should couple to any particular leptoquark. 

The formulae presented above we have implemented into the LESKO-YFS Monte 
Carlo program |^, We have performed the technical tests of the matrix elements 
as implemented in LESKO-YFS and compared them with analytical results on several 
approximated forms of the above matrix elements, e.g. for Z + 7 or only Z exchange ex- 
cluded. Agreement of 4-5 digits was always found. Later, QED corrections were extended 
(to our X-exchange, i.e. non t-channel Z, 7 interaction) according to the prescription ex- 
plained in Ref. |]19[ for the similar type of modification of the Monte Carlo KORALZ |^ 
in Ref. 
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We stress that the leptoquark-quark-lepton vertices which we have assumed in (|^), 
corresponding to the Ui,Ui, V2,V2, and f/3 spin 1 and R2, R2, Si, Si and 5*3 spin ex- 
amples in Ref. 10), are intended to be generic and not exhaustive: it is straightforward 



to include more general coupling scenarios into our LESKO-YFS calculational frame- 
work, should this become necessary. We point-out further that we may identify our 
states in models (1), ■ ■ ■ , (4) for spin 1,0 respectively with the corresponding charge and 
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mass eigenstates formed from the states Sl, T, Sr, Sr, Dl, Dr and D in the notation 



of Refs. |T^, |I6|; for, at scales ~ 200 GeV, we expect the SU2L x Ui EW symmetry 
to be broken with leptoquark states of the same charge and color mixed into the re- 
spective mass eigenstates and it is these mass eigenstates that we have used in (|1|). In 
other words, from the HERA data we we may have that eg resonates into the lepto- 
quark X, q = u, d, models (1) and (2) with t-channel X exchange in e~^e~ qq and 
F = in the language of Refs. [|12|, |15|, [l^, or that eq resonates into the leptoquark 
X, q = u, d, models (3) and (4) with u-channel X exchange in e^e~ qq and F = —2 
in the language of Refs. |]12|, |16[, where F is the fermion number of X. For definite- 
ness, we have assumed strong isospin symmetry for simplicity; it is trivial to relax this 
last assumption, should more data render this necessary. For completeness, we then 
note that the interaction Lagrangian densities which we used to arrive at the results (§) 

are, for S = 1, C^I^h^ = -i9x/V2) [xj^'^'^'^eHru'fj + xl-'^'^^^enrd^ + h.c.y and 

4tS!H^ = (xi-^/'^"47''e^, +xi"'/')"M-7^e^^, + /I.C.), where here H = L,R, 
H' = L when H = R and H' = R when H = L, and we have defined ipL(R) = Pl(r)'4^ 
as well as ip'^ = Cip'^ for all ip when T denotes transposition. The the chiral projections 
are here Pl = (1 — 7^)/2,Pr = (1 + 7^)/2 and C is the charge conjugation matrix in 
an appropriate representation. Thus, C^^t^^l^^ describes the interactions used to derive 
the result (Q) for Q = 5/3, S = 1, H = L, . . . , and so on. For S = 0, the analogous 
formulas follow from the corresponding replacements vector X x Dirac vector current 
<;=^ scalar X x Dirac scalar current. Finally, to facilitate contact with Refs. 16|, 
we note that, if we make the simplest possible assumption about the leptoquark weak 
isospin mixing matrix, that is that our spin 1 states are composed of only those states in 
Refs. |TB|, |T§, then we may identify {T^-^/3\ (T(-2/3) _^-(2/3)^^^| ^ ^^(.5/3)^ ^(-2/3) | 

in model (1), {^("^/s), ^-^Z^)} ^ {x(-5/3)^ x^-^/^)} in model (2), {d[-^/'\ D^''^''^} ^ 
{X(-4/3),x(-V3)} in model (3), and {d'^-^^'^\ D^-^/^^ ^ {X^-^/^), X^-^s)} in model 
(4), with the attendant coupling constant relations 2g = gx, \/2g = gx, V^g = gx, and 
= 9x, respectively where g is the coupling constant in Refs. |jl5|, [1^ (an analogous 
transformation holds for the spin case). 

In the Fig. 1, we exhibit the comparison of the expectations of the models (l)-(2) in 
comparison to the SM expectations for several observables available in the HERA data 



at low (2500 < < 15000 GeV^) for a value of 5 in the range allowed [|T3[ by 



the LEP2 data. Models (3) and (4) are simply too close to the SM expectation in this 
regime to describe the data for any value of 6 of order 1 in magnitude so that we do not 
bother to indicate these latter predictions explicitly in the comparisons with the data. 
Values of 6 several orders of magnitude in size larger than 1 are already excluded by the 
lack of observation of pronounced anomalies in the HERA e~ data for models (3) and 
(4) p. Shown are the and ye distributions in Fig. 1(a) and 1(b), respectively, for 
2500 <Ql< 15000 GeV^ and in Fig. 1(c) the distribution for 100 GeV < < 180 
GeV in the HERA kinematics in the hard scattering CM system for which we have the 
definitions s = {pp +pi)^, Q\ = -i, = Ql/{2{pi - P2)Pp), ye = Ql/{sXe), so that 
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Me = ^sXe is the invariant mass of the outgoing positron-jet system- the subscript e 
means that the observable is computed from variables on the e line |l|, |^. Here, is the 
incoming proton 4-momentum. Our e kinematics is defined so that all photons within a 
cone of half-angle 5° is included in the definition of the outgoing e four- vector, to give a 
more realistic simulation of the definition of the outgoing e in the HERA detectors. We 
also impose a cut on the initial state photon radiation by requiring that a fraction of 
the e energy lost to this radiation is < 0.1. We see that for gx such that 5 = —0.7 for 
example the = models are well within the experimental errors at HERA |^ of the 
SM expectations. By comparison, the S = \ models are disfavored at low y^,. 

We then show in Figs. 2 and 3 the comparison of the HI and ZEUS data, respectively, 
with the models for the observables Ue and Me for HI, and the observables xda and 
Qda (see Refs. [jl], || for their definition) for ZEUS, where for HI we show data for Ql with 
0.1 <ye< 0.9, for with Me > 180 GeV, for Me with Ql > 15000 GeV^, 0.1 < ye < 0.9, 
and for ye with Ql > 15000 GeV^ and for ZEUS we show data for xda with y^A > 0.25, 
Q'da > 5000 GeV^, and for Qfj^ in the regions of the anomalous numbers of events. 
The respective experimental acceptances [0, |^ are included in the above results. In all 
six distributions, we see that, for 6 = —0.7, there is reasonable agreement between the 
data and the theoretical prediction for models (1) and (2) for S = whereas for S* = 1 
the comparison with the data in Figs. 2a,b,d and 3a,b is less favorable but still tolerable 
(within 3a). Thus, we must await more data for a detailed check of the models (1) and (2). 

In our plots we have used the upper limit ^ value Tx = 2 GeV; if we use the theoretical 
lower limit Tx = 5'|^Mx/327r(487r) for S" = 0(1) the required value of S in Figs. 2, 3 changes 
to ~ —0.95 but our qualitative conclusions do not change. We have also checked that 



varying our structure functions between those in Ref. p2| and those in Refs. [21 1 does not 



affect our results in any significant way at the level of precision of the HERA data. What 
we can do is to emphasize that, if models (1) and (2) are indeed correct, we do not expect 
similar anomalies in the e~ data because that would involve scattering from the sea and 
this would be strongly suppressed. Further, we do not have to have charged current events: 
to get them, we need to add more terms to the interactions in the C\ni'^~H^ given above 
but in our general phenomenological framework, this is neither mandated nor forbidden. 
We also comment on the expectations at FNAL. What they are in general depends on 
the assumed BR to the initial state channel at HERA for the final state decay as well as 
on the detailed coupling scenario between gluons and X: as emphasized in Ref. [^], a BR 
near unity with only two higher multipoles in the gluon-X coupling scenario may already 
be disfavored by the DO+CDF [^, |2^, ^ limits if this combination can be done without 



unforseen systematics although as Ref. [18| has stressed the situation depends in detail on 



cuts, BR's, etc. and the implied signal may be small enough that it is missed and hence 
the matter is still unsettled; but, in our general phenomenological approach, we have no 
reason to expect a BR near unity or a naive coupling between leptoquarks and gluons: an 
infinite series of multipoles is in principle allowed in our effective low energy Lagrangian 
for X and gluons, with the truly underlying renormalizable theory awaiting yet higher 
energy probes to reveal itself. Thus, we can only encourage the FNAL experiments to 
continue their search for the observation of X as well. 
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Figure 1: Distributions of Mf. and ye obtained from LESKO-YFS for the Standard Model process 
(solid histograms) and in the presence of the scalar (full circles) and vector (open circles) leptoquarks 
in models (1) and (2) - both models give identical distributions. All the results were 
obtained for the MRSA parametrization of the proton structure functions. 

In summary, we have investigated the comparison with HERA data of the expectations 
on the leptoquark interpretation of the anomalous positron-jet phenomenon at DESY. 
We used the LESKO-YFS || Monte Carlo event generator so that higher radiative 
corrections to the Born level leptoquark signal are calculated with the YFS exponentiated 
0{a) LL /9o,i residuals, for both initial and final state radiation with finite-pT effects in this 
radiation properly realized, in the framework of Ref. [^, We find that, in agreement 
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Figure 2: Distributions of Ql, ye and Me obtained from LESKO-YFS for the Standard Model 
process (thin-line histograms) and in the presence of the scalar (thick-line histograms) and vector 
(dotted histograms) leptoquarks in models (1) and (2) with Mx = 200 GeV, Tx = 2GeV, 6s = 
—0.72, 6v = —0.68 (see the text for details) compared with the recent HI high events (★ symbols) 
as given in Ref. [ffl]. Av {v = Ql,ye, Me) denotes bin size for respective histograms. 



with the analysis in Ref. [0, a leptoquark couphng of < O.Sgw in our models (1) and (2), 
positron-quark resonances, where spin is preferred over spin 1, is in general agreement 
with the HERA data distributions at both high and low for Tx ^ 2GeV and that our 
models (3) and (4), positron-anti-quark resonances, like the SM itself, even in the presence 
of the rigorous higher order initial and final state radiation with finite-pT photon effects. 
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Figure 3: Distributions of xda and Q^/^ obtained from LESKO-YFS for the Standard Model 
process (thin-line histograms) and in the presence of the same scalar (thick-line histograms) and 
vector (dotted histograms) leptoquarks as in Fig. ^ compared with the recent ZEUS high events 
(o symbols) as given in Ref. H^v {v = xda,Q%a) denotes bin size for respective histograms. 

cannot account for the data for any reasonable value of the respective leptoquark coupling 
parameter. We look forward with excitement to more precise data. 
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